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ABSTRACT 


In  this  report  we  investigate  signal/noise  ratio 
enhancement  obtainable  on  single-station  LPE  recordings 
by  several  techniques.  The  processors  which  are  com- 
pared over  a data  base  of  20  events  include  simple 
band-pass  filtering,  match  filtering  with  a synthetic 
LR  signature,  and  time-varying  non-linear  processing 
which  examines  three- component  data  for  Rayleigh-wave 
characteristics.  For  the  suite  of  events  in  this 
study,  average  S/N  ratio  improvement  is  6 db  for 
match  filtering,  15  db  for  the  non-linear  PHILTRE 
processor  and  20  db  cascading  the  two  processors. 
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INTRODUCTION 


Recent  work  at  the  Seismic  Data  Laboratory  has 
shown  that  the  network  of  Long  Period  experimental  (LPL) 
stations  has  a 90%  threshold  of  approximately  4.7  (m^) 
for  detection  of  surface  waves  from  global  events 
(Von  Seggern,  1973) . This  is  not  nearly  low  enough  to 
insure  detection  of  surface  waves  from  many  events  of 
interest,  especially  underground  nuclear  explosions  of 
low  yield.  Since  these  stations  record  data  in  digital 
form,  it  is  relatively  easy  to  attempt  to  improve  this 
threshold  by  digital  processing.  Many  signal-enhance- 
ment processing  experiments  have  been  reported  for 
other  long-period  seismic  systems  (e.g.,  Capon  et  al.; 
1969;  Binder,  1970;  Simons,  1968);  but  because  of  the 
quite  different  response  of  LPL  systems  (compared,  for 
instance,  to  WWNSS  or  LRSM  lonb-period  responses)  we 
feel  that  a thorough  experiment  witii  available  pro- 
cessing methods  is  warranted.  In  this  report  we 
examine  and  compare  the  signal  enhancement  effected 
on  a suite  of  events  by  the  following  methods: 
band-pass  filtering,  match  filtering,  and  non-linear 
time-varying  processing. 
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METHODS 


Band-pass  Filtering 

As  a reference  point  for  judging  the  effectiveness 
of  our  inore  sophisticated  signal  enhancement  methods, 
a narrow  band-pass  filter  was  also  applied  to  the  LPli 
data.  We  employed  a program  (LOPAZ)  which  simulates 
analog  filtering  of  digital  data.  The  corner  frequencies 
were  set  to  .025  and  .040  Hz  (25  to  40  seconds  period) 
with  24  db  rolloff  on  each  side.  This  filter  should 
produce  an  output  which  will  be  free  of  any  microseismic 
noise  and  which  will  take  maximum  advantage  of  the  mini- 
mum in  earth  noise  at  30  to  40  seconds  period  (Savino 
et  al. , 1972) . 

Match  Filtering 

This  is  a signal  enhancement  technique  whereby  one 
crosscorrelates  a time  window  containing  the  desired 
signal  with  a reference  waveform  having  presumably  a 
nearly  identical  character.  Three  different  approaches 
are  commonly  used  in  constructing  this  reference  wave- 
form. One  is  an  empirical  approach  which  uses  a seismic 
signature  with  large  S/N  ratio  from  the  same  region  as 
the  signal  to  be  match- filtered  (Alexander  and 
Rabenstine,  1967;  von  Seggern,  1971).  A second  approach 
is  to  construct  a swept- f requency  or  chirp  waveform 
from  an  analytic  formula  (Capon  et  al.,  1969).  A 
third  is  to  synthesize  the  reference  waveform  from 
known  or  estimated  ~urface-wave  velocities  (Alexander 
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and  Rabenstine,  1967);  this  is  usually  done  by  inverse 
transforming  frequence-domain  data  points  defined  by 
the  phase  velocities  at  each  Fourier  frequency. 

For  this  study  a variation  of  the  third  approach 
was  employed.  Subroutine  SURFIL  (Kimball  and  Johnson, 
1970)  synthesized  the  reference  waveform  in  the  me 
domain  point  by  point.  It  uses  group  velocity  informa- 
tion provided  for  periods  between  50  and  10  seconds  by 
subroutine  SURTT  (Masse  and  Crowley,  1968;  and  Kimball, 
1969)  . This  in  global  grid  form,  each  grid  area  being 
assigned  one  of  sixteen  possible  Rayleigh-wave,  group 
velocity,  dispersion  curves;  SURTT  defines  the  great 
circle  path  between  the  epicenter  and  station  and  then 
averages  the  slownesses  at  each  period  over  the  multiple 
grid  areas  traversed.  The  envelope  of  the  waveform  con- 
structed by  SURFIL  is  shaped  by  the  instrument  ampli- 
tude responses;  however,  attenuation  was  not  taken 
into  account  in  constructing  the  waveform.  Adding  the 
attenuation  would  not  significantly  affect  results 
since  the  attenuation  factor  over  the  20  to  50  seconds 
band  is  not  highly  variable  (Tryggvason,  1965)  and  the 
precision  of  the  amplitude  spectrum  is  not  as  nearly 
as  important  for  match  filtering  as  that  of  the  phase 
spectrum  (Alexander  and  Rabenstine,  1967).  For  uni- 
formity a length  of  450  seconds  was  always  used  with 
this  synthetic  reference  waveform. 

The  program  MATCH  (Johnson,  1970)  which  calls  these 
subroutines  has  within  it  a band-pass  filter  which  windows 
the  frequency-domain  amplitudes  of  the  transformed  time 
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trace.  Its  corner  frequencies  are  .0165  and  .067  Hz; 
the  falloff  is  linear  to  zero  amplitude  at  .01  and  .10 
Hz,  respectively.  This  gives  an  effective  bandpass 
for  periods  between  approximately  12  and  80  seconds. 

This  band-pass  filter  was  used  in  all  the  match-filter 
processing  presented  in  this  report.  This  is  not  an 
effective  filter,  however,  for  most  of  the  LPE  record- 
ings used  here,  since  the  LPE  response  allows  little 
noise  energy  outside  this  band  to  appear  on  the 
recordings  anyway,  and  so  the  filter’s  contribution 
to  S/N  ratio  enhancement  in  the  match-filter  results 
in  probably  small  (2  or  3 db)  . 

Non-linear  Time-varying  Filtering 

Time-varying  adaptive  filters  using  three- component 
recordings  have  been  used  previously  for  seismic  signal 
enhancement.  Sax  and  Mims  (1965)  discussed  the  use  of 
polarization  filters  to  enhance  short-period  data.  Choy 
and  McCamy  (1973)  used  a polarization  filter  on  LPE  data 
recorded  at  OGD  to  inspect  the  relative  phase  angle 
between  vertical  and  horizontal  motion  in  order  to 
enhance  Rayleigh  waves.  They  also  used  separately  an 
azimuth  discriminator  to  enhance  Love  waves. 

For  this  study  the  time-varying  processor  used  is 
the  program  PHILTRE  developed  by  Simons  (1968)  . This 
program  operates  on  all  three  traces  after  rotating 
the  horizontal  traces  to  be  radial  and  transverse  to 
the  great-circle  path  from  epicenter  to  stations.  It 
is  thus  dependent  on  at  least  crude  epicenter  locations 
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for  the  event  being  processed  and  also  on  good  opera- 
tion and  calibration  of  all  three  instruments.  After 
the  traces  are  transformed  to  the  frequency  domain, 

PHILTRE  applies  three  criteria  to  pass  the  Rayleigh 
wave  signal  of  interest  and  attenuate  all  else.  First, 
the  frequency  components  are  weighted  according  to  the 
deviation  of  the  direction  of  approach  of  the  energy, 
as  given  by  the  relative  amplitudes  of  the  rotated 
horizontal  traces.  Second,  the  frequency  components 
are  weighted  according  to  whether  they  have  the  expected 
90°  phase  lag  between  vertical  and  radial  traces.  Third, 
they  are  weighted  according  to  whether  they  have  the 
expected  ellipticity  for  LR  motion  (values  of  0.6  to 
1.0  for  the  frequencies  of  interest).  A full  explana- 
tion of  the  scheme  can  be  found  in  Simons  (1968).  We 
used  weighting  exponents  of  8,  4,  and  0 for  the  azimuth, 
phase  lag,  and  ellipticity  criteria,  respectively.  ( ie 
ellipticity  criterion  was  suppressed  since  it  is  fre- 
quency and  station  epicenter  dependent,  and  we  did  not. 
have  specific  information  on  it.)  The  weighted  harmonic 
components  are  inverse  transformed  for  each  trace  to 
produce  the  output.  Using  data  sampled  at  0.5  pts/sec 
to  speed  up  processing  time,  256-second  time  windows 
are  transformed  to  the  frequency  domain.  These  time 
windows  are  successively  stepped  down  the  record  1/8 
of  the  window  length  (32  seconds)  to  produce  overlapping 
outputs  which  are  then  averaged  to  produce  a single 
continuous  output  for  each  component.  The  same  broad- 
band prefilter  as  used  in  the  match  filtering  is  use 
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in  the  PHILTRE  processing.  Although  gains  are  equa- 
lized between  the  three  components  for  computing  purpos 
they  are  not  equalized  on  the  plots  comprising  the  top 
six  traces  of  Figures  4 through  23. 

A Cascaded  Processor 

As  a final  processor  we  match-filtered  the  output 
from  PHILTRE:  that  is,  crosscorrelated  the  PHILTRE 

vertical  output  with  the  synthetic  waveform  used  pre- 
viously on  the  raw  vertical  trace. 


DATA 


For  this  study  the  receivers  are  at  three  of  the 
LPE  sites:  CTA  (Charters  Towers,  Australis),  KIP 

(Kipapa,  Hawaii)  and  KON  (Kongsberg,  Noway).  Twenty 
events  varying  from  3.7  to  5.5  (mean  = 4.61)  were 
randomly  chosen  witn  epicentral  distances  varying  from 
33°  to  137°.  Figures  1 through  3 show  location  of 
these  stations  and  their  associated  events;  and  coordi- 
nates and  magnitude  of  the  events  are  included  in 
Table  I.  Epicenters  are  from  the  NOS  List  with  the 
exception  of  3 from  the  LASA  daily  bulletin.  Most  of 
the  Rayleigh  waves  in  this  data  base  are  near  the 
visual  detection  threshold;  we  in  fact  know  that  two 
(KIP-8,  KIP-10)  were  undetected  in  routine  analysis  of 
the  LPE  data  at  the  SDL. 

The  frequency  response  of  the  LPE  system  is 
centered  at  approximately  40  seconds  period  to  minimize 
background  noise  and  hopefully  to  improve  detection  of 
earthquake  Rayleigh  waves  (Savino  et  al.,  1972). 

Seismic  data  are  digitally  recorded  at  the  sites  with 
a sampling  interval  of  one  second. 
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RESULTS 


Figures  4 through  23  are  30-minute  plots  for  each 
of  the  20  events.  Each  figure  shows  11  traces:  traces 

1-3  are  the  vertical,  radial  and  transverse  broad-band 
filtered  input  to  PHILTRE,  traces  4-6  are  PHILTRE  out- 
put, trace  7 is  the  synthetic  LR,  trace  8 is  the  cross- 
correlation of  traces  1 and  7,  trace  10  is  the  narrow- 
band  filter  output,  and  trace  11  is  the  raw  vertical 
component  shown  as  a reference.  Table  I lists  the  S/N 
improvement  in  db  of  the  various  processors.  Figure  24 
shows  histograms  of  the  enhancement  distribution  over 
the  20  events  using  increments  of  3 db.  S/N  ratio  is 
defined  in  this  report  as  a visual  measurement  of  the 
maximum  peak-to-peak  signal  amplitude,  divided  by  the 
rms  noise  amplitude  computed  from  a 256-second  window. 
The  noise  window  is  the  first  256  seconds  on  the  plots; 
its  position  is  held  invariant  for  each  input  or  output 
trace  for  that  event. 

Although  for  the  most  part  we  let  the  tabulated 
values  and  illustrated  traces  speak  for  themselves, 
a few  results  deserve  comment.  The  mean  processor  gain 
is  4 to  5 db  for  the  tight  band-pass  filter,  6 to  7 db 
for  match  filtering,  15  to  16  db  for  PHILTRE,  and  19 
to  20  db  for  PHILTRE  + MATCH.  Choy  and  McCamy  (1973) 
have  reported  to  a 6 db  average  gain  over  a unspecified 
number  of  cases  when  t\  ey  used  a Rayleigh-wave  enhancer 
operating  on  phase  angle  only.  Our  PHILTRE  results  are 
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much  better,  perhaps  reflecting  the  superiority  of 
Simons'  processor  and  the  fact  that  we  have  included 
in  the  PHILTRE  results  that  improvement  obtained  by 
band-passing  the  data  with  a broad-band  filter.  Table  I 
lists  three  events  for  which  the  improvement  of  the 
cascaded  processor  is  less  than  for  PHILTRE  alone  and 
four  events  where  the  cascaded  processor  does  better 
than  the  sum  of  the  two  individual  processor  results. 
These  results  are  not  so  surprising  when  we  consider 
we  are  really  trying  to  relate  results  of  MATCH  on 
PHILTRE  output  to  results  of  MATCH  on  raw  data,  the 
discrepencies  arising  from  the  time-varying  aspects 
of  the  non-linear  PHILTRE  processor. 

Another  observation  is  that  events  which  are 
greatly  enhanced  by  one  processor  were  seldom  so 
greatly  enhanced  by  another.  For  example,  the  events 
match  filtering  did  best  on  had  only  average  S/N 
enhancement  in  PHILTRE,  while  most  of  the  events  with 
high  S/N  improvement  in  PHILTRE  (over  20  db)  showed 
little  match-filter  improvement  over  that  of  the  tight 
band-pass  filter.  One  is  led  to  suspect  that  this  is 
due  to  the  differing  characteristics  of  the  processors. 
Match  filtering  assumes  that  all  group  velocity  grid 
values  are  fairly  well  known  and  neglects  multipath 
effects  (Binder,  1970).  PHILTRE  is  unaffected  by 
multipathing  from  back  azimuths  near  the  great  circle 
azimuths,  but  is  hampered  by  extraneous  coherent  energy 
approaching  from  the  same  azimuth  as  the  signal  to  be 
enhanced. 
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Considering  enhancement  as  a function  of  travel 
path,  the  majority  of  PHILTRE  gains  over  20  db  are 
attributable  to  circum- Pacific  events  arriving  at  KIP 
while  those  with  hign  S/N  improvement  in  match  filtering 
include  only  some  of  these,  and  two  from  the  North 
Atlantic  to  KIP  across  the  United  States.  Events  with 
poorer  than  average  S/N  improvement  on  both  PHILTRE 
and  MATCH  include  events  whose  travel  paths  were  from 
Peru  to  CTA,  Greenland  to  KIP,  and  the  North  Atlantic 
to  KO.'I. 

It  is  important  to  know  whether  a processor  will 
do  as  well  on  weak  signals  ac  on  strong  signals.  To 
check  this,  the  processors’  S/N  improvement  populations 
were  divided  into  "high- gainer"  and  "low-gainer" 
classes  to  see  if  such  a division  correlates  with  in- 
put S/N  ratio.  The  histograms  of  S/N  ratios  of  the 
raw  data  shown  in  Figure  25  demonstrates  that  there 
are  at  least  as  many  "high-gainers"  among  the  weak 
signals  as  among  the  strong  signals  (in  fact  both 
MATCH  and  PHILTRE  appear  to  do  better  on  weak  signals 
for  this  population).  For  comparison,  a visual  detec- 
tion threshold  of  1:1  (S  /N  ) is  shown  relative  to 

PP  PP 

the  event  distributions.  We  do  not  have  signals  well 
below  this  threshold  and  therefore  cannot  determine 
how  rapidly  the  processors  degrade  with  low  S/N  ratios. 
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CONCLUSIONS 


Over  our  suite  of  20  events,,  match  filtering 
improved  S/N  ratio  by  6 to  7 db . This  figure  is  roughly 
"'hat  is  expected  (Capon  et  al.,  1969);  and  its  2 db 
improvement  over  band-pass  filtering  alone  on  these 
events  is  arso  roughly  what  has  been  achieved  before 
(Binder,  1970).  Although  match  filtering  shows  some 
2 db  advantage  over  band-pass  filtering  here,  PHILTRE 
shows  a very  significant  11  db  advantage  over  it. 

I he  combination  of  PHILTRE  + MATCH  achieved  on 
the  average  a remarkable  S/N  ratio  improvement  of  20 
db  over  raw  traces  (this  includes  the  wide  band-pass 
prefilter  as  described  above).  We  can  comfortably 
predict,  though,  that  this  20  db  cannot  hold  up  for 
signals  below  the  visual  threshold.  The  PHILTRE 
processor  depends  upon  finding  some  signal  component 
with  which  to  reject  the  noise;  and  when  noise  begins 
to  dominate,  the  processor  may  never  find  the  signal 
component.  The  calculations  required  by  PHILTRE  + 

MATCH  are  not  cumbersome,  and  real-time  implementation 
would  be  easily  attainable.  However,  disadvantages  of 
these  processors  are  that  application  of  PHILTRE 
requires  all  three  components  and  knowledge  of  back 
aximuths  to  events  being  processed  and  that  applica- 
tion of  MATCH  requires  at  least  crude  estimates  of  the 
epicenter  locations  of  events  being  processed.  These 
event  location  parameters  should  be  readily  obtainable 
from  short-period  network  detections  though. 
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SUGGESTIONS  FOR  FURTHER  STUDIES 


Three  further  studies  are  suggested  by  the  results 
of  this  report:  one  would  be  a comparison  of  signal- 

to-noise  ratio  obtainable  for  an  LPE  site  versus  an 
array,  using,  for  instance,  KON  and  NORSAR  data  for 
common  events.  The  cascaded  PHILTRE-MATCH  processor 
would  operate  on  the  KON  data  while  the  NORSAR  data 
would  first  be  beamformed  before  the  same  cascaded 
processor  was  applied  to  it.  A second  study  would 
use  the  cascaded  processor  to  see  how  far  the  detec- 
tion threshold  of  an  LPE  station  can  be  lowered  by 
‘ processing  the  data  at  predicted  Rayleigh-wave  arrival 

times  for  visually  undetected  events  of  a given  low 
mb  for  a specific  region.  A third  study  would  be  to 
evaluate  what  capability  PHILTRE  has  to  separate 


mixed  events  recorded  by  a single  three-component  site 
due  to  its  azimuthal  discrimination  and  phase  isola- 
jp’  tion  characteristics. 
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Figure  5.  Processing  Results  for  CTA- 


Figure  7.  Processing  Results  for  CTA- 
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Figure  9.  Processing  Results  for  CTA- 
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Figure  12.  Processing  Results  for  KIP- 
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Figure  15.  Processing  Results  for  KIP- 
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Figure  17.  Processing  Results  for  KIP- 


Figure  18.  Processing  Results  for  KIP- 


Figure  19.  Processing  Results  for  KIP- 10 
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Figure  22.  Processing  Results  for  KON- 


Figure  23.  Processing  Results  for  KON- 
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Figure  24.  Histograms  of  signal  enhancement  by  each  of  the  four 
processors  on  25  events 
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figure  25.  Distribution  of  raw  signal/noise  ratios 
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